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Abstract 

We present some theoretical results relevant to the direct dark matter detection experiments, paying par- 
ticular attention to directional experiments, i.e. experiments in which, not only the energy but the direction of 
the recoiling nucleus is observed. In directional experiments the detection rate depends on the angle between 
the line observation and the sun's direction of motion. Since, however, the direction of observation is fixed with 
respect the earth, while the Earth is rotating around its axis, in a directional experiment the angle between the 
direction of observation and the Sun's direction of motion will change during the day. So the observed signal in 
such experiments will exhibit a very interesting and characteristic periodic diurnal variation. 
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> ■ 1 Introduction 

m The combined MAXIMA-1 Q], BOOMERANG 2 , DASI [3] and COBE/DMR Cosmic Microwave Background 



(CMB) observations [4 imply that the Universe is flat [5] and that most of the matter in the Universe is Dark [BJ, 
i.e. exotic. These results have been confirmed and improved by the recent WMAP data [7]. Combining the the 
data of these quite precise experiments one finds: 

fl b = 0.0456 ±0.0015, n CD M = 0.228 ±0.013, fl A = 0.726 ± 0.015. 



Since any "invisible" non exotic component cannot possibly exceed 40% of the above Vic dm [Ej, exotic (non 
\T2 \ baryonic) matter is required and there is room for cold dark matter candidates or WIMPs (Weakly Interacting 
j_j ■ Massive Particles). 

Since the WIMP's are expected to be very massive and extremely non relativistic, they are not likely to excite 
the nucleus. So they can be directly detected mainly via the recoiling of a nucleus (A,Z) in elastic scattering. In the 
standard nuclear recoil experiments, first proposed more than 30 years ago [S], one has to face the problem that the 
reaction of interest does not have a characteristic feature to distinguish it from the background. So for the expected 
low counting rates the background is a formidable problem. Some special features of the WIMP-nuclear interaction 
can be exploited to reduce the background problems. Such are: i) the modulation effect: this yields a periodic 
signal due to the motion of the earth around the sun. ii) backward-forward asymmetry expected in directional 
experiments, i.e. experiments in which the direction of the recoiling nucleus is also observed, iii) transitions to very 
low excited states, if present in the target: in this case one need not measure nuclear recoils, but the de-excitation 
7 rays, (iv) detection of electrons produced during the WIMP-nucleus collision or subsequent X-rays. 

In our previous work [lOj [TT] [12] we have found that the observed directional rate is characterized by a strong 
dependence on the angle between the direction of nuclear recoil and the direction of the sun's motion, with respect 
to the center of the galaxy. The apparatus will, of course, be oriented in a direction defined in the local frame, i.e. 
the line of observation will point to a point in the sky specified, in the equatorial system, by right ascension a and 
inclination 5. Thus due to the rotation of the earth in a directional experiment the angle between the direction of 
observation and the Sun's direction of motion will change during the day. So, since the event rates sensitively depend 
on this angle, as we have recently suggested |13j one may be able to observe a diurnal variation. In the present 
paper we will explore this novel feature of directional experiments, i.e. the variation of the data with a period of 
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24 hours. This variation is independent of the Earth's rotational velocity. Its amplitude depends, in addition to 
some of the usual parameters specifying the standard WIMP event rate, on the inclination of observation, i.e. the 
angle between the direction of observation and the axis of the Earth's rotation. Those features cannot be masked 
by any known background. 



2 Standard (non Directional) Rates 

The event rate for a given WIMP velocity v is given by [T3] 



dR P * mt -V^^A^ n ( itV dt ^' v) , (1) 



du m x Am p V Mp / du 

where m t / Am p is the number of nuclei in a target of mass m tl A is the mass number of the target, u is the energy 
transfer in some convenient units, u — Q /Qq , Qq = 40 A~ 4 / 3 MeV , and 

^ = -^=^2^' (2) 
du V< v 2 > 2{^ r bvy 

with F(u) the nuclear form factor. We must now fold dt(u, v)/du with the velocity distribution f(y, £) in the local 
frame. Here we use y — v/vq, v$ is the sun's velocity with respect to the center of the galaxy, and £ is the cosine 
of the the WIMP velocity with the polar axis (chosen to be the sun's direction of motion) . We will assume in this 
work that the velocity distribution is Maxwell-Boltzmann in the galactic frame, namely 

/(y') = ^e- (y ' )2 , y' = -, (3) 

where v' is the WIMP velocity in the galactic frame. Obviously to obtain the distribution in the local frame one 
must take into account the motion of the Sun and Earth: 

y' = y + z + S (sin ax — cos a cos jy + cos a sin jz) , (4) 

where 6 — Vi/vo = 0.136 with v\ the velocity of the earth around the sun, a the phase of the earth (a=0 on June 
3rd) and 7 ~ w/6. 

In the case where 6 = 0, we get for the differential (with respect to the energy transfer) rate 
dR p x m t ( Hr\ 2 1 9 — /i 2 J~2 



-a n (^j 2 V7^A 2 ^ ^- = ^a 2 F 2 ( U )^ (aV^), (5) 
*o(a^)= [ V " C ydy 2n [ d£/(y,0, (6) 



du m x Am 



where 

a = (V2n r bv )- 1 . (7) 

We note that the function ^>q(x), a decreasing function of the energy transfer, depends only on the velocity 
distribution and is independent of nuclear physics. An additional suppression as the energy transfer increases 
comes, of course, from the nuclear form factor F(u). 

Note that the nuclear dependence of the differential rate comes mainly from the form factor and, to some extent, 
via the parameter a. Integrating the differential event rate from u m i n = Eth/Qo, which depends on the energy 
threshold, to u max = yl sc /a 2 (where y esc = 2.84 since v esc = 2.84uq) we obtain the total rate 

R=^^a n (^) 2 V^v^A% t=r a *du^. (8) 
m x Am p VMp/ J Umin du 

In the case where 5 ^ we proceed as above and finally get for the differential rate an expression similar to ([5]) 
where now one has to replace: 

tfofc) -> *o(a0 +*i(a;,7,'5)cosa, (9) 

where 4 , i(x, 7, S) is obtained as the linear term in an expansion of the velocity distribution in powers od 5. Thus 
one has to replace: 

dt dr dt dh 

> — = — + — cosa, (10) 

du du du du 
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m x — > GeV 

Figure 1: The quantity h, entering in the coherent mode discussed in this work, is shown as a function of the 
WIMP mass in GcV for Q mm = for a light system, 32 S (left) and 127 I or 131 Xe (right). 
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Figure 2: The coherent total event rate, as a function of the WIMP mass in GeV, is exhibited assuming a nucleon 
cross section a n = 10~ 43 cm 2 =10 _7 pb for Q m in — for a light system, 32 S (left) and a heavy system, 131 Xe (right). 



, , i r max , dft 

t — > til + h cos a), a = — j du — . 

t ./,, . du 



(11) 



Thus we finally get 



R = 



Px m t ( Mr 



m x Arrip \fi p 



— V< v 2 >A 2 t(l + h cobol) 



(12) 



3 Directional Rates 

In this instance the experiments will attempt to measure not only the energy, but the direction of the recoiling 
nucleus as well [HI [TTJ [121 E] • Let us indicate the direction of observation by: 

e = (e x ,e y ,e z ) — (sin cos sin 8 sin cos 8). (13) 

Now Eq. ^ must be motified since no integration over the azimuthal angle specifying the direction of the outgoing 
nucleus is performed. Thus Eqs. ([1]) and ([2] ) now become 



du Jdtr m x Am p \fi p J 2tt V du , dir 

F ( u )^—[-T? 6 [v-e- —7— . (15) 



dw / dir V< v 2 > 2(n r bv) 2 V f/fo, 

The factor 1 / (27r) enters since we are going to make use of the same nuclear cross section as in the non directional 
case. 

The total event rate is obtained after folding with the velocity distribution and integrating the above expression 
over the energy transfer. It can be cast in the form: 



R dlr = fx^^f^Y^^^K^^ 

m x Arrip \ n p J Zir 



l + /i c (8,$)cosa + /i s (8,$)sina;) . (16) 
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Figure 3: Due to the diurnal motion of the Earth the angle 8, between the line of observation and the sun's 
direction of motion, is changing as the earth rotates. We show the angle 6 as a function of a, essentially time, for 
various inclinations 5. The intermediate thickness, the short dash, the long dash, the fine line, the long-short dash, 
the short-long-short dash and the thick line correspond to inclination 6 = — 7r/2, — 37r/10, — 7r/10, 0, 7r/10, 37r/10 
and 7r/2 respectively. We see that, for negative inclinations, the angle O can take values near it, i.e. opposite to 
the direction of the sun's velocity, where the rate attains its maximum. 



We found it convenient to use a new relative parameter 43 as well as the parameters h m and ctQ given by: 

rc(6,$) = — , h m (Q,$) cos (a + a ) = /i c (6, $) cosa + /i s (6, $) sina. (17) 

These parameters depend, of course, on the reduced mass /z r , the WIMP velocity distribution and, to some 
extent, on the nuclear physics via the nuclear form factor. The parameter n gives the retardation factor of 
the directional rate, over and above the factor of l/(27r), compared to the standard rate. Since, however, the 
unmodulated amplitude is independent of <fr one can integrate over $ so that the suppression of 1/ (27r) drops out 
for this term. Because of the existence of both cos a and sin a terms the time dependence will be of the form 
cos (a + ao), with the phase ao being direction dependent. Thus the time of the maximum and minimum will 
depend on the direction. In other words the seasonal dependence will depend on the direction of observation. So 
it cannot be masked by irrelevant seasonal effects. 

Before concluding this section we like to consider the case connected with the partly directional experiments, 
i.e. experiments which can determine the line along which the nucleus is recoiling, but not the sense of direction 
on it. The results in this case can be obtained by summing up the events in both directions, i.e. those specified 
by (8, $) as well as (tt — 9, $ + it). A given line of observation is now specified by 0, $ in the range: 

< 6 < tt/2, < $ < tt. 

The parameter h entering the non directional case is shown in Fig. [T] Using the above parameters t and 
employing Eq. (|12p we obtain the time independent total event rate. The obtained results are shown in Fig. [5] 
The total modulation, since it is defined relative to the time independent part, is still given by Fig. [TJ 

4 Diurnal Variation 

Up to now we have considered the event rate in a directional experiment in fixed direction with respect to the 
galaxy in the galactic system discussed above. The apparatus, of course, will be oriented in a direction specified in 
the local frame, e.g. by a point in the sky specified, in the equatorial system, by right ascension a and inclination!! 
8. This will lead to a diurnal variation of the event rate [T51 114) . 

The galactic frame, in the so called J2000 system, is defined by the galactic pole with ascension <Si = 
12' 1 51™ 26.282 s and inclination Si = +27° 1 42.0l" and the galactic center at a 2 = I7 h 45 m 37.224 s , 
$2 = —(28° 56 10.23 ). Thus the galactic unit vector y, specified by (oti,Si), and the unit vector s, specified 

1 This is done to avoid or minimize, e.g., the dependence on the parameters of the particle model and, in particular, the unknown 
WIMP-nucleon cross section. Furthermore the directional experiments can also obtain the rate in all directions. Thus k maybe of 
experimental interest. 

2 We have chosen to adopt the notation a and 6 instead of the standard notation a and <5 employed by the astronomers to avoid 
possible confusion stemming from the fact that a has already been used to designate the phase of the Earth and 8 has been used for 
the ratio of the rotational velocity of the Earth around the Sun by the velocity of the sun around the center of the galaxy 
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Figure 4: The time dependence (in units of the Earth's rotation period) of the parameter n for various in- 
clinations 5 in the case of the target CS2. The intermediate thickness, the short dash, the long dash, the 
fine line, the long-short dash, the short-long-short dash and the thick line correspond to inclination S — 
— 7r/2, — 37t/10, — 7t/10, 0, 7r/10, 37r/10 and ir/2 respectively. Due to the diurnal motion of the Earth different angles 
O are sampled as the earth rotates. 
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Figure 5: The same as in Fig. @]in the case of the Iodine target. The intermediate thickness, the short dash, the 
long dash, the fine line, the long-short dash, the short-long-short dash and the thick line correspond to inclination 
5 = -tt/2, -371-/10, -tt/IO, 0, tt/10, 3tt/10 and tt/2 respectively. 



by (cb,^), can be expressed in terms of the celestial unit vectors i (beginning of measuring the right ascension), 
k (the axis of the Earth's rotation) and j = k X i. One finds 

y = -0.868i - 0.198? + 0.456A; (galactic axis) , 

x = —s = 0.055z + 0.873j + 0.483A; (radially out towards the sun) , 

z = xxy = 0.494* - 0.445j + 0.747fc (the sun's direction of motion). (18) 

Note in our system the x-axis is opposite to the s-axis used by the astronomers. Thus a vector oriented by (a, 5) 
in the laboratory is given in the galactic frame by a unit vector with components: 

/ y \ / -0.868 cos a cos 5 -0.198 sin & cos <5 + 0.456 sin 6 \ 

I x = 0.055 cos a cos 5 + 0.873 sin a cos £ + 0.4831 sin 5 , (19) 
\ z J \ 0.494 cos a cos 6 - 0.445 sin a cos 5 + 0.747 sin 6 / 

where ao = 282.25° is the right ascension of the equinox, 7 w 33° was given above and 6p = 62.6° is the angle the 
Earth's north pole forms with the axis of the galaxy. 

Due to the Earth's rotation the unit vector (x, y, z), with a suitable choice of the initial time, a — ao = 2ir{t/T), 
is changing as a function of time 

The angle O scanned by the direction of observation is shown, for various inclinations <5, in Fig. 3. We see that 
for negative inclinations, the angle can take values near n, i.e. opposite to the direction of the sun's velocity, 
where the rate attains its maximum. 

The equipment scans different parts of the galactic sky, i.e. observes different angles 0. So the rate will change 
with time depending on whether the sense of the recoiling nucleus can be determined along the line of recoil. 
The results depend, of course, on the WIMP mass and the target employed. We will consider a light (the time 
dependence of k is exhibited in Fig. [4]) and an intermediate-heavy target (the time dependence of k is exhibited in 
Fig. E). 



5 Discussion 



In directional experiments one measures not only the energy of the recoiling nucleus but also the direction the 
nucleus is recoiling. Some of the requirements that should be met by such detectors have recently been discussed 
[TBI [17] . To fully exploit the advantages of such detectors one should be able to distinguish between recoils with 
momenta p and — p (sense of direction), which now appears to be feasible. Some of the predicted interesting features 
of directional event rates persist, however, even if it turns out that the sense of motion of recoils along their line of 
motion cannot be measured. Such experiments, given a sufficient number of events, provide an excellent signature 
to discriminate against background. One expects, of course, a smaller rate by observing in a given direction. In 
the most favored direction, opposite to the sun's direction of motion, the event rate is rs down from that of the 
standard non directional experiments, if a specific angle $ is chosen. Since, however, n is independent of the angle 
<£, one can integrate over all azimuthal angles and thus the retardation is only of order k (k ~ 1 in the most favored 
direction). The parameter k is essentially independent of any particle model parameters other than the WIMP mass. 
It does depend, however, on the assumed velocity distribution and to the nuclear form factor. Finally we have 
seen that in directional experiments the relative event rate for detecting WIMPs within our galaxy, as given by the 
parameter n, will show a periodic diurnal variation due to the rotation of the Earth. The time variation is larger 
in the case of light WIMP and/or light target. So from this perspective the lighter target is preferred. The time 
variation arising after the inclusion of the modulation parameters Qi m ,a) or (h c ,h s ) is expected to be even more 
complicated, since these parameters depend on both and One expects a diurnal variation on top of the annual 
variation characteristic of the usual modulation effect entering both directional and non directional experiments. 
Such effects will be discussed elsewhere. 
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